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Abstract 


In this work, we consider a collocation method for solving the pantograph-type Volterra Hammerstein integral equations 
based on the first kind Chebyshev polynomials. We use the Lagrange interpolating polynomial to approximate the 
solution. The convergence of the presented method has been analyzed by over estimating for error. Finally, some 
illustrative examples are given to test the accuracy of the method. The presented method is compared with the Legendre 
Tau method. 
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1 | Introduction 


QO) This paper aims to obtain the exact approximations by collocation method based on the first kind 
Li , . : srg 

cons Chebyshev polynomials for solutions of the following pantograph Volterra Hammerstein integral 
Computational ¥ 

Algorithms and equations, 


Numerical Dimensions. 
This article is an open qs 
access article distributed | z(s) = g(s) + f R(t, s)G(t, z(t))dt, 0<q<1, se [01]. (1) 
under the terms and 0 


conditions of the Creative | Where g(s), R(s, t) and G(t, z(t)) are smooth functions on their domains. We assume that the function 


Commons Attribution ; 5 : k oF ; š 
" G(t, z(£)) satisfies in the Lipschitz condition with respect to the second variable. 


(CC BY) license 
(http://creativecommons. , i . . 
org/licenses/by/4.0). These equations appear in many branches of science such as control theory, biology, ecology and etc 


[1]-[7]. Severral methods have been proposed to solve the integral equations in [8]-[11], [5], [12]-[16], 
[6], [17]-[26]. Ansari and Mokhtary [27] presented the Legendre Tau method for solving Eg. (1) and 
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discussed applying the spectral methods to obtain a reliable numerical solution for Eg. (1) according to 
the well-known existence and uniqueness theorems in [27]. 


Here, we approximate the solution of Eq. (7) using the Lagrange interpolating polynomial. We consider 
a collocation method based on the first kind Chebyshev polynomials for solving Eq. (1) and get a non- 
linear system that can be solved by Newton method to obtain the solutions at the grid pionts. In order 
to use the theory of orthogonal polynomials, we transform the Eg. (7) to [-1,1] and analyze the 
convergence of the presented method. Some numerical examples prepared to test the efficiency and 
accuracy of the proposed method. We compare the numerical results of the presented method with the 
Legendre Tau method in [27]. 


2 | Algorithm of the Method 


The first kind of Chebyshev polynomials Ty(s) are orthogonal at [-1,1] with respect to the weight 
function w(s) = (1 - s*)-¥? and are defined by [28]: 


Tn(s) = cos(Ncos “1(s)). 

All Chebyshev polynomials Ty(s) can be generated by the following recursion relation 
Tn(s) = 2sTy-1(8) = Tyn-2(8), N= 2,3, ssi 

With 


To(s) = 1, Ti(s) = 8S. 


The Gauss quadrature formula 


N 
f wads > £(5,)We- 
| k=0 


is exact for any polynomial of degree < 2N + 1. 


Now, in order to use the theory of orthogonal Chebyshev polynomials, applying the following change 
of variable 


t=dned, -1<n<Ç 0<q<1, 


a 1<ÇC<1 
s = — =; -1<C<l. 
2 2 

Por Eq. (1), we have 


Gs cae 
20) = BO + | RCWC(n,2))an, -1< 6 <1. 2 


~ 1 1 
RGM =3R(50+5, 


Using the Lagrange interpolating polynomial, we can approximate Z(C) as 
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IZO) = Şi 1,(Q2(s;). (3) 


j=0 
Where 
1m) = T= , i= O(N. 
j=0,j+i 
(2j+1)r 


N 
And s; = -cos ( )} are the Chebyshev Gauss quadrature points. Also, from Eg. (3), we can write 
j=0 


2N+2 


G(1S(2@)) = $i (OG (s; z(s;)). a 
Now, collocating the nodes s;, i= 0(1)N, in Eq. (2), we have 
Si — ~ 
2(s;) = &(s;) + f R (sn) G (n, 1§(2)) dn, i=O()N. (5) 
-1 
Solving the non-linear system (5) by the Newton method, we get the values of 2(s j)- Substituting these values 


into Eq. (3), the values of REO) ate obtained, for all € [-1,1] . 


3 | Convergence Analysis 


In this section, using some useful inequalities and lemmas from [29], [30], we will try to determine over 


estimate for |II§(2(2) - ZOllrzc1,1) 


Theorem 1. Assume that REO) is the spectral collocation solution of Fg. (2) given by Eq. (3). If the 
functions (C), R, n) and G(n,2(n)) ate sufficiently smooth and G(n,2(n)) satisfies in the Lipschitz 


condition with respect to the second variable and e(C) = REO) — 2(C), then we have 


lelg can SCN™ (Bhan, yy + Lenan) 


a © 
+ CNF (|G || +N-* LZ, 24 n) 
Proof. The Eq. (3) can be written as follows 
Si — ~ Si — ~ 
us) =86)+ | Rem G(n 182) an+ f Rum Gln, zm) 4 
g 7 (7) 


-f R (sın) G(n, 2(n)) dn. 
Multiplying Eg. (7) by l;(Č) and sum up from 0 to N, we have 
INZ) = 1f(&©) 
+ IĘ( f R (Gn) (z (n 1S%)) - G(n,20))an) 8 
+IÑ f R n) G(n,2(n)) dn). 
Subtracting Eq. (8) from the Eg. (2), we have 
e(C) = Lz + 1S i RG) (z (n,1S(2m)) - Ela,20)} An) + Ls: ) 


Where 


L = 18(&) - 80, 
Cas = Gr eo 
Ls SINJ RCo G(n,2(n)) an- | RG) G(n,2(n)) d 
The Eq. (9) can be rewritten as follows 
delat a R (Gn) (z (n,1S()) - (nz) dn +L; + L4. (10) 


Where 


18 f N (n1S(m)) - E(n,20n))a n) 
- [Rew (Efike) - Ehr z)]an 


To simplicity, in the following the notation ||. || has been used instead of |]. Ilz2, EEI Now, from the Eq. 
(10), we have 


lle(Oll < IL2ll + Il f Ren (E (n 1S%))- G 20) ) nel + al (11) 


Now, in order to determine over estimate for error, we should find the over bound for right hand side 
of the inequality Fg. (17) term by term. 


Using inequality (5.5.22) from [29], we have 


IL2ll < CN |B enan (12) 

Where 
m 1/2 
2 =~ (KINI 
bean =f Dy e (13) 
k=min(m,N+1) 
Por Eg. (12), using Eg. (13) with m =1 and Lemma 3.8 from [30], we have 
: a(R) 
ILI] < CN | f S G(nZn))dn + RG OG(G ZO) | 
z (Gn) 14 

< CN | f aR) H ) G(n, zanl + IRC OME o) a 

< CN*(C]G(n, zw) + IRE ONCE 20). 
Thus, from Eq. (14), we have 

IILgll < CN-1||GI]. (15) 
To find over bound for ||L4|, similar to process of obtaining Eg. (15), we have 

Lull < CN“ [|G (n 1S(2p)) - G(n,20)]}. (16 


Since G(n, 2(n)) satisfies in the Lipschitz condition with respect to the second variable. Thereby from 
Eq. (16), we have 
|[Lall < CNLllell. (17) 


Where L is the Lipschitz constant. Now, applying (5.5.22) from [29] for Eg. (17), we have 
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Using Lipschitz condition and (5.5.22) from [29], we have 
Ce = Ciz | ae Sas 
IS Rn) (z (n 1S()) = E(n,200)) dn || < CLN Zhan 1)° 
Combining the above estimates with Eq. (11) leads to the intended error estimate Eq. (6). 
4 | Numerical Results 
In this section, we use the proposed method for solving some numerical examples to test the accuracy of 
the method. We get the numerical results by Wolfram Mathematica 12.2. The presented method has been 


compared with the Legendre Tau method in [27]. For simplicity, we use the notation ||. || instead of |]. llz2- 


Example 1. Consider the following pantograph-type Volterra Hammerstein integral equation [27] 


z(s) = g(s) + [TRE s)G(t,z(t))dt, s € [0,1]. 


Where 
1 
g(s) =- ge P(s)(s — sin(s)) + tan(s), 
-i7 
= 9710" 


R(s, t) = exp(s) cos?(t), 
G(t,z(t)) = z(t). 


And the exact solution is z(s) = tan(s). The function G(t,z(#)) satisfies in the Lipschitz condition with 
respect to the second variable. 


Table 1 contains a comparison of the L?-norms of the error function between the presented method and 
the Legendre Tau method in [27] for Example 1 with q = > Figs. 1 and 2 show the plots of the obtained 
errors |z — IX(z)| by the presented method in Example 1 for N = 8 and N = 32, respectively, with q = z, Fig. 
9 shows tending ||z — ISI to zero by increasing N in Example 1 with q = > 
Table 2 contains the obtained errors ||z — ISI by the presented method for Example 1 with q = = Figs. 3 
and 4 show the plots of the obtained errors |z — I§(z)| by the presented method in Example 1 for N = 8 and 
N = 32, respectively, with q = Z, Fig. 10 shows tending ||z — IS(z)|| to zero by increasing N in Example 1 


. 7 
with q = 5; 


Example 2. Consider the following pantograph-type Volterra Hammerstein integral equation [27] 
qs 
2(s) = g(s) + | Rit-s)G(t2()dt, s € [0,1], 
0 
Where 
1 3 3 
g(s) = 3 (4 cosh (=) — 4 cosh(s) + s sinh E) + sinh È), 
ai 
q= 4 
R(s, t) = cosh(s — t), 
G(t,z(t)) = sinh(z(t)). 


And the exact solution is z(s) = sinh À). The function G(t,z(#)) satisfies in the Lipschitz condition with 


respect to the second variable. 


Table 3 contains a comparison of the L?-norms of the error function between the presented method and 
the Legendre Tau method in [27] for Example 2. Figs. 5 and 6 show the plots of the obtained errors 
|z — IS (z)| by the presented method for Example 2 with N = 8 and N = 32, respectively. Fig. 17 shows 
tending ||z — ISIl to zero by increasing N in Example 2. 


Example 3. Consider the following pantograph-type Volterra Hammerstein integral equation 
qs 
(3) = g(s) + | Rit-s)G(t2()dt, s € [0,1], 
0 
Where 


2 4 2 
g(s) = 2.125 — 0.42s? + cos È) — 2 cos (=) — 0.125 cos (=) —1.8s sin (=) = 


. {As 
0.225s sin (=). 
5 
8 
eT 


R(s,t)=s+t, 

G(t,z(t)) = z(t)*. 
And the exact solution is z(s) = cos(=). The function G(t,z(t)) satisfies in the Lipschitz condition with 
respect to the second variable. 


Table 4 contains the obtained errors ||z — ISIl by the presented method for Example 3. Figs. 7 and 8 
show the plots of the obtained errors |z - I§(2)| by the presented method in Example 3 for N = 8 and 
N = 32, respectively. Fig. 12 shows tending ||z — ISI to zero by increasing N in Example 3. 


Table 1. Comparison of the presented method and the Legendre Tau method in [27] for example 1 
1 


with q= 7 
N jiz- IS(z)l| forthe |z- IS(ll for the 
Presented Method Legendre Tau Method 

2  4.99e -02 == 

4 1.97e -03 6.17e — 02 
& 1:19e=05 7.96e — 03 
16 2.03e -10 1.62e — 04 
32 1.296 =15 8.66e — 08 
64  8.66e -16 3.31e — 14 


Table 2. The obtained errors ||z — I§(z)|| by the presented method for example 1 with q = Z 


N 2 4 8 16 32 64 
liz- IS(@Il 6.01e -02 1.50e-03 1.23e-05 1.97e-10 877e= 16. 1.74e— 16 


Table 3. Comparison of the presented method and the Legendre Tau method in [27] for example 2. 


N |z- IS(z)| for the Presented Method _||z — IS (Il for the Legendre Tau Method 


2  4.23e -13 = 

4  1.69e -14 6.29e — 04 
8  8.49e -16 1.17e — 05 
16 5.69e -16 1.24e — 08 


32 3.95e -16 4.94e — 14 
64 1.04e —16 = 
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Table 4. The obtained errors ||z — I$ (z)|| by the presented method for example 3. 


N 2 4 8 16 32 64 
lz- IȘ(z)|| 2.85e -04 1.29e-07 1.35e—13 8.30e-14  1.62e—-15 1.01e -15 


Error 


8.x 10 16t 
6.x10 1®---------- 
4.x10 16} 
2.x10 1% 


02 04 06 08 107 
Fig. 2. Plot of the obtained errors |z — I§(z)| by the presented method for N = 32 in example 1 with q = > 
Error 
8.x10°} 
6.x10 ê 
4.x10 °] 


2x10% 


Fig. 3. Plot of the obtained errors |z — IÇ(z)| by the presented method for N = 8 in example 1 with q = = 


Error 


Fig. 5. Plot of the obtained errors |z — I§(z)| by the presented method for N = 8 for 


Sj 


02 04 06 08 1.0 


example 2. 


Fig. 6. Plot of the obtained errors |z — I§(z)| by the presented method for N = 32 for example 2. 


Sj 
0.2 04 06 08 10 


Error 


1x10 18h -------------- 
8.x 10° 14 
6.x10'4--------------- 
4.x10 14 
2.x10 14} 


Fig. 8. Plot of the obtained errors |z — I§(z)| by the presented method for N = 32 for example 3. 
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Fig. 9. Graph of tending ||z — IĘ(z)|| to zero by increasing N in example 1 with q 
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Fig. 12. Graph of tending ||z — IĘ(z)|| to zero by increasing N in example 3. 


5 | Conclusion 


In this work, we applied a spectral collocation method to the numerical solution of the pantograph-type 
Volterra Hammerstein integral equations based on the first kind Chebyshev polynomials. We employed 
the Lagrange interpolating polynomial to approximate the solution. 


The convergence of the presented method is analyzed. Some numerical examples are prepared to test 
the accuracy of the proposed method. We observed that the proposed method is more accurate than 
Legendre Tau method in [27]. The other polynomials can be used for solving these equations such as 
Chelyshkov polynomials, Legendre polynomials, other kinds of Chebyshev polynomials, and so on. 
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